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ABSTRACT 



ConXexX. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the 
mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. 
/K\ms. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks 
around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. 

Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distribu- 
tions and to predict the spectral slopes at (sub-)mm wavelengths. 

Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our 
models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only 
be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This 
dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes 
would become too low) or due to efficient conversion of dust into larger, unseen bodies. 

Key words, accretion, accretion disks - circumstellar matter - stars: formation, pre-main-sequence - infrared: stars 
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1. Introduction 

Circumstellar disks play a fundamental role in the formation 
of stars as most of the stellar material is believed to be trans- 
ported through the disk before being accreted onto the star 
(Lynden-Bell & Pringle 1974). At the same time circumstellar 
disks are thought to be the birth places of planets. Understanding 
the physics of circumstellar disks is therefore the key to some of 
the most active fields of astrophysical research today. 

However, observing these disks in order to learn about the 
physical processes taking place in their interior is a challenging 
task. Strom et al. (1989) and Beckwith et al. (1990) were the first 
to use observations at mm-wavelengths to confirm that many of 
the observed pre-main sequence (PMS) stars showed excess ra- 
diation above the spectrum of a T Tauri star. While these single 
dish observations provided valuable insight in dust masses (since 
mm observations probe not only the thin surface layers, but the 
bulk of the dust mass in the mid-plane), recent developments in 
the field of (sub-)mm interferometry grant the possibility to con- 
strain models of disk structure and evolution of protoplanetary 
disks by fitting parametric models to the observed radial profiles 
(e.g., Andrews et al. 2009; Isella et al. 2009). Spatially resolving 
the disks is important since it ensures that low millimeter spec- 
tral slopes are not just an artifact of high optical depth. 

Today, mm spectral slopes are known for quite a num- 
ber of disks and spatially resolved observations indicate that 
the low values measured in these samples are related to grain 
growth (e.g., Testi et al. 2003; Natta et al. 2004; Rodmann et al. 



2006). Grains are believed to collide and stick together by van 
der Waals forces, thus forming larger and larger aggregates 
(Dominik & Tielens 1997; Poppe et al. 2000; Blum & Wurm 
2008). Due to this loose binding, collisions with velocities in 
excess of a few m s"' may lead to fragmentation of the aggre- 
gates. 

Larger samples of radially resolved mm spectral slopes 
are expected in the near future, but still no study so far in- 
terpreted mm observations using simulated grain size distri- 
butions but rather used simple parametric power-law distribu- 
tions with an upper size cut-off. In this work, we use a state 
of the art dust grain evolution code (similar to Brauer et al. 
2008a; Birnstiel et al. 2010) to derive steady-state grain distri- 
butions where grain growth and fragmentation effects balance 
each other. We self-consistently solve for the grain size distribu- 
tions and the disk structure to predict fluxes at mm wavelengths 
and the radial dependence of the mm spectral index. Comparing 
these results to observed values in the Taurus and Ophiuchus 
star-forming regions allows us to test predictions of the theory 
of grain growth/fragmentation and to infer constraints on grain 
properties such as the critical collision velocity and the distribu- 
tion of fragments produced in collision events. 

Grains orbiting at the Keplerian velocity in a laminar gas 
disk feel a constant head wind (caused by the gas rotating 
slightly sub-keplerian) which forces them to spiral inwards 
(Weidenschilling 1977). If this drag is as efficient as laminar 
theory predicts, all dust particles which are necessary to explain 
the observed spectral indices, would quickly be removed (see 
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Table 1. Parameters of the model grid: Mdisk is the total disk 
mass, Q-t is the turbulence parameter, Uf is the critical collision 
velocity, /vac is the grain volume fraction of vacuum and ^ is the 
index of the distribution of fragments (see Eq. 3). The parame- 
ters of the fiducial model are highlighted in bold face. 
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Brauer et al. 2007). In this Letter, we therefore assume that ra- 
dial drift is halted by an unknown mechanism. Under this as- 
sumption, we find that low values of the mm spectral index can 
be explained by the theory. We show that in order to explain the 
observed flux levels, the amount of observable' dust needs to 
be reduced by either reducing the dust-to-gas ratio (perhaps by 
radial drift at a intermediate efficiency or during an earlier evo- 
lutionary epoch) or by dust particle growth beyond centimeter 
sizes. 



2. Model description 

2.1. Disk model 

We consider disks around a PMS star with a mass of 0.5 Mq, 
bolometric luminosity of 0.9 Lq and effective temperature of 
4000 K, at a distance of 140 pc, which are typical values for 
the sample of low-mass PMS stars studied in Ricci et al. (2010, 
hereafter RIO). To derive the disk structure we adopted a mod- 
ified version of the two-layer models of passively irradiated 
flared disks developed by Dullemond et al. (2001) (following 
the schematization by Chiang & Goldreich 1997), in which we 
have relaxed the common assumption that dust grain proper- 
ties are constant throughout the disk. For the disk surface den- 
sity we adopted the self-similar solution for a viscous disk 
(see Lynden-Bell & Pringle 1974) with parameters lying in the 
ranges observationally constrained by Andrews et al. (2009). 
The surface density gradient y and the characteristic radius 
(for the definitions, see Hartmann et al. 1998) are assumed to be 
J - I and /?c = 60 AU, respectively. Throughout this work we 
assume a constant dust-to-gas mass ratio of 1%. 

2.2. Dust model 

We use a coagulation/fragmentation code as described in 
Brauer et al. (2008a) and Birnstiel et al. (2010) to simulate the 
growth of dust particles. Particles grow due to mutual col- 
lisions (induced by Brownian motion and by turbulence, see 
Ormel & Cuzzi 2007) and subsequent sticking by van der Waals 
forces. We assume the dust particles to be spheres of internal 
density ps and vary ps to account for porosity effects. However, 
we do not treat a dynamic porosity model (see Ormel et al. 2007; 
Zsom & Dullemond 2008). 

With increasing collision velocity Am, the probability of 
sticking decreases and fragmentation events start to become im- 



' by "visible" or "observable" dust we mean the dust particles which 
are responsible for most of the thermal continuum emission at (sub-)mm 
wavelengths which are typically smaller than a few centimeter in radius. 



portant. In this Letter, we use the fragmentation probability 

( if Am < Mf - 6u 

Pf^ll if Am > Mf (1) 



where Mf is the collision velocity above which particles are as- 
sumed to fragment and 6u is the transition width between coag- 
ulation and fragmentation (taken to be 0.2 Mf). Recent studies of 
collision experiments (Guttler et al. 2010) and numerical simu- 
lations (Zsom et al. 2010) indicate that there is also a regime in 
which particles may bounce. However, since this topic is still not 
well understood, we will omit these effects in this work. 

Radial drift is an, as yet, unsolved problem (Brauer et al. 
2007; Birnstiel et al. 2009). However there are several ef- 
fects such as spiral wave structure (e.g., Cossins et al. 2009), 
density sinks (e.g., Brauer etal. 2008b) or zonal flows (e.g., 
Johansen et al. 2009) which may reduce the effectiveness of ra- 
dial drift. In this paper we assume that radial drift is ineffective 
since we focus on the question whether observations can be ex- 
plained through the physics of grain growth and fragmentation. 
The question to answer in this case is not, how to retain parti- 
cles at these radii, but rather how to create them there in the first 
place. 

To investigate this problem, we simulate the physics of par- 
ticle growth and fragmentation until a steady state between both 
processes develops. Since the relative velocities for particles typ- 
ically increase with grain radius, we can relate the fragmentation 
velocity to a certain grain size (which defines the "fragmentation 
barrier", see Birnstiel et al. 2009) 

2Sg M^ 

^max ~ ' (2) 

above which particles fragment (with Eg, at and being the 
gas surface density, the turbulence parameter and the sound 
speed, respectively). Mf and are assumed to be radially con- 
stant with at values within a range expected from theoretical (see 
Johansen & Klahr 2005; Dzyurkevich et al., 2010, in press) and 
observational works (see Andrews et al. 2009). Grains which 
reach flmax will experience high velocity collisions causing them 
to be eroded or even completely fragmented. The resulting frag- 
ments can again contribute to growth processes at smaller sizes 
and the grain size distribution will at some point reach a steady 
state where gain and loss terms caused by coagulation and by 
fragmentation cancel out at all sizes. 

Particles will need a certain time to grow to the fragmenta- 
tion barrier. The time to reach the steady state will therefore be a 
few of these growth time scales. Depending on the distance to the 
central star, the steady state is typically reached after a few thou- 
sand years at 1 AU up to about 1 Myr at 100 AU. The mean ages 
of the sources in our sample are » 2 Myr and a; 0.5 - 1 Myr for 
the Taurus and Ophiuchus PMS stars, respectively. Since radii 
around 40 - 80 AU dominate the observed emission at (sub-)mm 
wavelengths, we expect most of the samples to be in or at least 
close to a steady state. 

If the highest collision velocity that turbulent motion induces 
(depending on at and Cj) is smaller than the critical collision ve- 
locity Mf, then (at least some) particles do not fragment (i.e. the 
break through the fragmentation barrier) and a steady state is 
never reached. Due to this scenario, some of the possible com- 
binations of the parameter values (see Table 1) do not reach a 
steady state and are therefore not included in the results. 

The shape of the steady state grain size distributions is influ- 
enced mainly by five parameters: the previously mentioned at. 
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Mf, £g, the temperature T (through the sound speed Cs), and by 
the prescription of fragmentation. In our models, we assume the 
distribution of fragments to follow a power-law number density 
distribution, 

n{m) oc rrT^, (3) 

with an upper end at nif. We consider fragmentation and crater- 
ing, as described in Birnstiel et al. (2010). Recent experiments 
suggest^ values between 1.07 and 1.37 (see Guttler et al. 2010). 
In this work, we consider ^ values between 1 .0 and 1 .8. 

To calculate the dust opacity of a given grain size distribu- 
tion we adopted the same dust grain model taken in RIO, i.e. 
porous composite spherical grains made of astronomical sili- 
cates, carbonaceous materials and water ices (see RIO for the ref- 
erences to the optical constants). The ratio between the fractional 
abundance of each species comes from Semenov et al. 2003 and 
models with three different porosities have been considered in 
this Letter (see Table 1). We used the Bruggeman mixing the- 
ory to combine the refractive indices of the different materials 
and to calculate the dust opacity of the composite grains. The 
opacity induces probably the largest uncertainties in our calcu- 
lations as grain composition, grain structure and temperature ef- 
fects may lead to largely different opacities (see, for example 
Henning & Stognienko 1996). 

2.3. Comparison to observations 

We compare the (sub-)mm SED generated by our models with 
observational data of RIO and Ricci et al. (in prep), more specif- 
ically the flux at 1 mm (/^imm) and the spectral index between 
1 and 3 mm iF{A) oc A^"'-^'™). The samples considered in- 
clude all the class II disks in the Taurus-Auriga and p-Oph star- 
forming regions respectively for which both the central PMS 
star and the disk are observationally well characterized through 
optical-NIR spectroscopy/photometry and (sub-)mm photome- 
try/interferometry. To calculate the dust opacity as a function of 
wavelengths and radius, and the temperature in the disk mid- 
plane, we iterated the two-layer disk model (keeping the profile 
of 2g constant in time) with the dust model described above until 
convergence is reached. Once the physical structure of the disk 
is determined, the two-layer disk models return the disk SED 
which can be compared to the observations. 

The influence of the different parameters on the calculated 
Q'i-3mm values can mostly be understood by a simple model for 
a dust distribution, as used in RIO (cf. Fig. 3 in RIO): for max- 
imum grain sizes much smaller than the observed wavelengths, 
the spectral index of the dust opacity j8i-3mm {k{A) oc A^^"'^'") is 
constant, while it decreases for flmax values larger than sub-mm 
sizes. In between (at a few tenth of a mm), there is a peak which 
is caused by an increased opacity of grains with sizes similar 
to the observed wavelength. The relation between a'i-3mm and 
/?i-3mm depends on the emitting spectrum and the optical depth. 
For a completely optically thin disk in the Rayleigh- Jeans regime 
j6i-3mm - a'i-3mm - 2. Howcvcr, if the emitted spectrum deviates 
from the Rayleigh-Jeans limit, then ySi-snun ^ ai-3mm - 2. In our 
models, Q'i-3mm -jSi-3mm tums out to be typically between 1 .4 and 
1 .7 if flniax is outside of the peak of opacity. 

3. Results 

3.1. Sub-mm fluxes and spectral indices 

For all possible combinations of the parameters shown in 
Table 1, we solved for the steady-state grain size distributions 
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Fig. 1. Influence of the parameters at (top left), fragmentation 
velocity (top right), disk mass (bottom left) and grain porosity 
(bottom right) on the observed fluxes and spectral indices. The 
black circle denotes the fiducial model whose parameters are 
given in Table 1 . The grey area represents the region in which 
the observed sources lie (see Fig. 2). 



and derived the ai-3n,m and Fimm values. As noted before, some 
of the models do not result in a steady state and are, therefore, 
not shown here. 

The top left panel of Fig. 1 shows the influence of the turbu- 
lence parameter a^. According to Eq. 2, the maximum grain size 
increases if at decreases. Depending on where a^n^ lies with re- 
spect to the opacity peak (see RIO, Fig. 3), Q'^3mm can increase 
or decrease with increasing at. In the simulations presented here, 
flmax is typically so large that increaseing at predicts larger spec- 
tral slopes. 

flmax is more sensitive to Uf (cf. top right panel in Fig. 1): 
the maximum grain size a^ax is proportional to m^, therefore a 
change of Uf by a factor of about 3 significantly changes ai-3,nm 
by increasing the grain size by about one order of magnitude. 
However many models with a fragmentation velocity of 10 m/s 
never reach a steady state. It is therefore not possible to explain 
lower Q'i-3mm valucs by a further increase of t/f alone. 

The way that Mdisk influences Fimm and ai-3mm is twofold. 
Firstly, a decrease in Mdisk (assuming a constant dust-to-gas ra- 
tio and a fixed shape of the disk surface density i.e. not varying 
Rc and y, see Section 2.1), reduces the amount of emitting dust 
and thus Fi,nm- Secondly, such a reduction in gas mass also re- 
duces flniax (Eq. 2), which tends to increase aumm- This com- 
bined trend is seen in Fig. 1. Hence, in order to explain faint 
sources with low ai-3mm, the amount of emitting dust has to be 
reduced while the disk gas mass stays large. This effect could 
be achieved in two ways; the amount of dust could be reduced 
by radial drift at a reduced rate (full radial drift would quickly 
remove all mm-sized grains, see Brauer et al. 2007) or only the 
"visible" amount of dust is reduced if some of the dust is al- 
ready contained in larger bodies. This latter case is predicted by 
our non-steady-state distribution models and will be discussed 
in more details in a forthcoming paper. 

In general, lower values of ^ translate to shallower grain-size 
distribution, which result in lower values of y6i-3mm (see Draine 
2006). The lower right panel in Fig. 1 does not seem to indicate 
a strong dependence on ^, however lower values of ^ (around 1 ) 
seem to be closer to the observations especially at large fluxes. 

Figure 2 shows the areas which are covered by our sets of 
simulations for different porosities in comparison to the obser- 
vational samples. It can be seen that only the brightest sources 
are covered by the simulations. The trend of larger Qri-3mm for 
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Flmm (mJy) 



Fig. 2. Observed fluxes at mm-wavelengths of the Taurus (black 
dots) and the Ophiuchus (blue dots) star-forming regions (see 
RIO and Ricci et al. 2010, in prep) and the areas covered by the 
simulation results for different vacuum fractions of the grains 
(varying all other parameters according to Table 1). 

larger vacuum fraction seems to be in contradiction with Eq. 2, 
since smaller grain volume density leads to larger Omax- However 
in this case, the opacity is much more affected by changing the 
grain structure: reducing the grains vacuum fraction increases 
the spectral index at sub-mm while it is reduced for longer wave- 
length. Therefore opacity effects outweigh the smaller changes 
in flmax- A more thorough analysis of opacity effects is beyond 
the scope of this letter, however it seems implausible that the 
large spread in the observations can be explained by different 
kinds of grains alone (see Draine 2006). 

3.2. Radial profiles of the dust opacity index 

The presented models also compute ffi-smm as function of ra- 
dius. From the point of view of comparison to observations, this 
is somewhat premature since observational methods are not yet 
able to provide reliable radial profiles of ai-amm (e-g-, Isella et al. 
2010, Banzatti et al., in prep.). However, the predicted radial 
dependence of y6i-3mm (shown in Fig. 3) agrees with the obser- 
vations so far It can be seen that the shape of most models 
looks similar, sUghtly increasing from j6i_3mn,-values around 0.5 
at 10 AU up to around 1.5 at 100 AU. The reason for this is 
that flmax depends on the ratio of surface density over tempera- 
ture. Under typical assumptions, a^ax will decrease with radius. 
An upper grain size which is decreasing with radius and stays 
outside the peak in the opacity results in y6i-3mm increasing with 
radius (cf. Fig. 3 in RIO). If the radially decreasing upper grain 
size fljnax reaches sizes just below mm, then the peak in opacity 
will produce also a peak in the radial profile of ySi-smm (the size 
of which depends much on the assumed opacity), which can be 
seen in Fig. 3. Thus, even though fln,ax is monotone in radius, 
jSi-3mm does not need to be monotone. 

4. Discussion and Conclusions 

In this Letter, we present the first in-depth comparison of simu- 
lated grain size distributions and observed mm spectral indices 
of YSOs in the Taurus and the Ophiuchus star-forming regions. 
Additionally we present the first predictions of the radial profile 
of the dust opacity index at mm wavelength which are consistent 
with the limits set by Isella et al. (2010). 
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Fig. 3. Predicted profiles of the dust opacity index at mm- 
wavelengths for different variations of the fiducial model. The 
colors correspond to the parameters shown in Fig. 1. 

Low values of the observed mm-slopes are quite naturally re- 
produced by our models, favoring low values of ^ and at as well 
as fragmentation threshold velocities above 1 m s However, 
a simultaneous match to the observed range of flux levels re- 
quires a reduction of the dust mass by a factor of a few up to 
about 30. This over-prediction of fluxes cannot be fixed by sim- 
ply reducing the disk mass since the predicted ai-3mm would be 
too large for smaller disk masses. Opacities induce a large un- 
certainty in the flux levels. However, considering the results of 
Draine (2006), it seems implausible that the large spread in ob- 
served fluxes for different disks with similar Q;i-3mm (which is 
probably even larger as very faint disks are not contained in the 
sample) can be explained by different grain mineralogy alone. 

The aforementioned reduction of observable dust could be 
due to radial drift at a reduced rate or during an earlier epoch 
(drift has been artificially suppressed in this work in order to ex- 
plain the low values of ai-3nim by > 1 mm sized grains). Another 
possible explanation is grain growth to even larger sizes as these 
bodies have a small opacity coefficient per unit mass. 

Finally, a different dependence between Q'i-3mm and the ob- 
served flux Fimm might also originate from disk surface densities 
profiles that differ from what we have assumed in this work. This 
possibility, as well as the effect of a different dust composition, 
will be considered in a future work. 
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